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Abstract 

We extend the picture of a transfer of nuclear spin-i polarization along a homoge- 
neous one-dimensional chain with the XY-Hamiltonian to the inhomogeneous chain with 
alternating nearest neighbour couplings and alternating Larmor frequencies. To this end, 
we calculate exactly the spectrum of the spin-^ XY-Hamiltonian of the alternating chain 
with an odd number of sites. The exact spectrum of the XY-Hamiltonian is also applied 
to study the multiple quantum (MQ) NMR dynamics of the alternating spin-^ chain. MQ 
NMR spectra are shown to have the MQ coherences of zero and ± second orders just as 
in the case of a homogeneous chain. The intensities of the MQ coherences are calculated. 

PACS numbers: 05.30.-d, 76.20.+q 

1 Introduction 

The discovery of the exact solution of spin-^ homogeneous one-dimensional chains with the 
XY-Hamiltonian PP , [2] gives the observable and measurable features unraveling the NMR 
dynamics of spin-^ homogeneous one-dimensional chains [3]. Although most of the NMR 
experiments for one-dimensional spin chains appear to be well-explained by means of the 
nuclear spin dynamics on the homogeneous spin chains, NMR spin dynamics beyond the 
homogeneous chains has attracted the attention recently. For example, the experiments 
4 demonstrate the propagation of spin wave excitations along the inhomogeneous spin 
chains, the mesoscopic echo has been observed [5] due to the reflections of the spin waves 
at the boundaries of the chain. 

The presented paper is aimed at exploring the key differences of the NMR of the 
homogeneous spin-i chain from the inhomogeneous spin-i chain. It may be difficult, 
however, to calculate the NMR responses for an inhomogeneous spin chain with a random 
variation of the nearest neighbour (NN) dipolar coupling. Thus, as a first step to unravel 
of the inhomogeneous effects, we treat the one-dimensional chain with an alternating 
spin-^ NN dipolar coupling and alternating Larmor frequencies. 

In the approach pQ, [2], the basic tool in exploring the ID spin-^ XY-Hamiltonian is 
the Jordan- Wigner transformation of the original spin-i XY-Hamiltonian to the Hamil- 
tonian of the free fermions. In this way, the exact thermodynamics for alternating infinite 
chains with spin-| the XY-Hamiltonian was explored in [BJ. The report [7j presents the 
exact spectrum of the XY-Hamiltonian with alternating couplings on the finite rings. 
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However, to our knowledge, the spectrum of the XY-Hamiltonian with alternating cou- 
plings on the open chains is lacking by now. This is what will be addressed in section 2. 
The derived exact spectrum of the spin-^ XY-Hamiltonian of the alternating open chains 
permits one to explain the transfer of the nuclear polarization along the alternating chains 
in section 3 and to calculate the MQ intensities of alternating spin-^ chains in section 
4. The concluding section 5 draws the distinction of the NMR dynamics on alternating 
chains from that on the homogeneous ones. 

2 Exact spectrum of spin-^ XY-Hamiltonian with 
alternating couplings on open chains 
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Figure 1: An open chain of the odd number, N, of spins with the alternating NN coupling constants 
D\ and D2 and alternating Larmor frequencies uj\, u>2- 

In this section, we derive the exact spectrum of the spin-^ XY-Hamiltonian 

N N-l 

H = ^ U n I nz + ^ D ntH+ x (ln,xln+l,x + In,yln+l,y) (1) 
n=l n=l 

of the open chain of the odd number, N, of sites with the alternating NN coupling con- 
stants D% and T)% and alternating Larmor frequencies u\ and UJ2, see Fig. 1. 

The nuclear spins are specified by the spin-| operators I na at sites n = 1, . . . , N with 
projections a = x,y, z. The Jordan- Wigner transformation 

n=n-\ \ 



In,— — In,x Hn,y — ( 2) j I[ z J C n , 

V 1=1 J 

In,+ = I n ,x+iln,y = (-2) n - 1 ( [J ^4, 



1=1 



2 



c n c n 2 ' 



(2) 



from spin-^ operators I na to the creation (annihilation) operators c+ (c n ) of the spinless 
fermions takes the Hamiltonian (^Q) into the Hamiltonian 



N 



i\ i 



U n C n C n 



n=l 



or in the matrix notations as 



N-l 



+ 2 E D n,n+1 {cnCn+1 + C+ +l C n } , 



(3) 



n=l 



N 



H = -c+ (D + 20) ^ - - J2 "n- 



(4) 



n=l 



In Eq. Q, we denote the row vector c + = (cf , . . . , c^), the column vector if = 
(ci, . . . ,Cjv)* (the superscript t represents the transpose) and specify the matrices f2 and 
D as 



(5) 
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Diagonalization of the matrix D + 2f2 is performed by the unitary transformation 

D + 2n = UAU + , A = diag{Ai,...,Ajv}, 
so that the new fermion operators 7^ and 7^ introduced by the relations 

N N 
n=l n=l 

bring the Hamiltonian (jlj) into the Hamiltonian 

N 1 N 



(6) 



(7) 



1 1 



(8) 



k=l 



71=1 



with energies of the free fermion waves. 

To go further in explicit calculations it is necessary to find the eigenvalues X u and 
eigenvectors \u v ) = (ui u , u 2v , ■ ■ ■ ,unu) 1 of the matrix D + 2ft, 

(D + 2Sl)\u v ) = \ v \u v ). (9) 

In Eq. ©, the components u n) „ at the even sites n = 2, 4, . . . , N — 1 obey the equations 

N-l 



Diu 2k -i, u + 2u 2 u 2k:U + D 2 U2k+l,v = Ku 2 k,u, k = 1, 2, 



(10) 
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We now fix the indices k and v in Eq. (|1U[) and eliminate from Eq. H10|) the components 
U2k-i,u and U2k+i,u at the odd sites (2k — 1) and (2/c + 1) by writing down the equations 
for the U2k-i,v, u 2 k+i,v from Eq. @ as 

D2U2k-2,u + 2u>lU2k~l,jy + D\1l2k,v = ^vU2k~l,v 

+ D 1 U 2 k+2,u = KU2k+l,w (11) 

Substituting into (jlUj) . we get the relations for the amplitudes u n ^ at the even sites, 

gU2k-2,u + gU2k+2,u = U2k,u, k = 1,2, . . . , (N — l)/2, (12) 

with the spatially independent coupling constant 

= Dl£2 n „v 

5 (A„-2u/i)(A„-2w2)-.D?-.Df I J 

The derivation of (|12() . <|13|) is a characteristic way in the theory of the real space renor- 
malization [S| which eliminates the half of the degrees of freedom belonging to the odd 
sites yielding the field equations on the (even) lattice sites with doubled lattice constant. 
Eq-s. (|T2*)) . (|T3|) govern the amplitudes U2 iV , U4,v, ■ ■ ■ , un-i,v ( including the amplitudes 
u 2,v, u N-i,v at the border sites 2 and N — 1 ) if we introduce the additional sites n = 
and n = N + 1 and put there 

u 0iU = 0, u N+ i )U = 0. (14) 

Eq. ()14j) implies the cutting off the lattice at the sites 1 and N, thus, preventing the 
fermions from escaping the lattice shown on Fig. 1. The solution of Eq. ()12|) conditioned 
by Eq. (JT1J) read 

U2k,u = A u sm[——), k,u = 1,2,..., (N-l)/2, (15) 



N + l 



with eigenvalues 



= wi+w 2 ± V(wi -£j 2 ) 2 + £>? A u , v = l,2, 



2-KV \ 2 



N- 1 



A, = l + 2,5 cos (j^- J + 5 = D 2 /D X . (16) 

Eq. (|16() gives 2^-^- = N — 1 eigenvalues A;, since for each index ^ = 1, 2, . . . , JV and for 

each superscript (+) and (— ) the relationship a£ = X^]i_ v holds. It is convenient to 
arrange (N — 1) distinct eigenvalues X u as 



A = i Ai +) , forz/ = l,2,...,(iV-l)/2 

fari/ = (JV + 3)/2,(JV + 5)/2,...,JV 



In accord with the enumeration of the eigenvalues A^ in Eq. (|17|) . the missed A^-th 
eigenvalue \ u stands for the index v = + 1 = (recall that N is odd). To find 
the eigenvalue Ajv+i and the eigenvector u jv+i at even sites n, the use is made of the 

2 2 

properties (see Eq. (jT^ ) 

un+±=0, n = 0,2,4,... ,N -1,N + 1. (18) 

n, 2 
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By the condition u N _ l n+\ = 0, the iV-th equation from Eq. @ gives immediately the 
sought eigenvalue 

A iv+i = 2u x . (19) 

2 

It remains to find the amplitudes u, n+i at the odd sites. The sought amplitudes obey 

k, 2 

the relations 

N + 1 , N - 1 

•Dl«2fc-l,i> + D 2 U 2 k+l,v = 0, i/= — - — , fc = l,2, — - — . (20) 

Eq. (|2*U1) gives the components of the eigenvector u jv+i belonging to odd sites n (up to 
the normalization coefficient B), 

u n> N±i =B-{-8)^, n = l,3,...,N. (21) 

Calculating the normalization coefficients A v in Eq. ()15|) and the coefficient £? in Eq. 
(|21|) . one finds all iV eigenvectors ju^) and all iV eigenvalues \ v of the Hamiltonian 



wi + W2 + V (wi - ^2) + f 1 i/ = l,2, 



jV-1 
2 



A„ = { 2wi, ^ = ^ • ( 22 ) 

Wl+W2 _y / ( Wl _ u;2 )2 + ^ A„ I/= ^,JV±5,...,JV 

For all indices ^ = 1,...,N except the index v = ^p-, the eigenvector \u u ) has the 
elements 

Au^[Ssm(^) +S in(^)], j = 1, 3, 5, . . . , N 
A v sn{0), j = 2,4,...,iV-l 

and the normalization coefficient 

A = 2K~2^| 1 

By Eq. (|2*T|) . the elements of the eigenvector \u n+i ) read 

2 

Vm = { 5 - ( - 5)V ' J = ?' 3 ' 5 '--;i iV . • (25) 



2 



0, j = 2,4,...,iV-l 



/ «5 2 -l \3 



with the normalization coefficient 

3 Spin-wave propagation in open chains with al- 
ternating couplings 

The spectrum of Eq.-s l|22[) - (j26|) of the XY Hamiltonian (JHJ) can now be applied to 
describe how the polarization at a single site of the alternating chain with equal Larmor 
frequencies uj n changes over the time. 
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Figure 2: Time course of the polarization (h,z)(t)/{h,z)(0) of the first spin, Eq. ({30]) . for the chain 
of 201 spins interacting with coupling coefficients D\ = 2tt ■ 4444 s _1 and D2 = 2ir ■ 6666 s" 1 and equal 
Larmor frequencies. The initial polarization is on the first spin. Insert shows the time course of the 
polarization on an earlier time interval < t < 4 ms. 



Let the initial polarization is on the single site j, hence, the spin dependent part of 
the initial density matrix at the high temperature approximation is, see PJ, 



p(0) = I jjZ = cfcj - -. 



Given the initial density matrix Q27|) and the Hamiltonian H 
mann equation (h = 1) 



.dp 



(27) 

the Liouville-von Neu- 

(28) 



is solved in terms of the fermion operators as follows, 
Pit) 



+ u j,i u j,m,e 2 



ll 7m- 



(29) 



Kl,m<N 



Denoting the polarization at j'-th. spin at the time moment t by (Ij' tZ )(t) and invoking 
Eq. I|27|). we get 



(Ij',z)(t) _ tr{p(t)I jl>z } 
<^>(0) " trWj 



J2 u luUf,ue * 



A,// 



Kv<N 



(30) 



Fig. 2 shows the time course of the polarization {Ii z )(t)/{Ii,z)(fy at site f = 1 when the 
initial polarization is also at site 1 and D2 = f-Di- Just as in the case of the polarization 
dynamics on the homogeneous chain with equal couplings, D2 = D\, [3], the dynamics 
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of the polarization of the alternating chain can be regarded as the propagation of the 
spin wave packet starting at site 1 and bouncing back and forth at the chain ends To 
calculate the return time £r for the wave packet to reappear at site 1, we, first, determine 
the group velocity of the waves described by the Hamiltonian H (jHJ) with the dispersion 
law T^Aj, (flB^I and equal Larmor frequencies, uj\ = u>2- By specifying the wave vector 
p = 2ttu/(N + 1), v = 1, 2, . . . , (N — l)/2, the dispersion law of Eq. (|T6*|) written down as 
+ 2<5cos(p) + 5 2 allows one to calculate the sought group velocity 



v = max <, a 

p 



d\(p) \ 
2 dp] 



(31) 



where a is the lattice constant and the coupling constant D\ is the minimal coupling 
constant among the two coupling constants D\ and Z?2- Thus, for N = 201-chain with 
coupling constants D\ = 27r-4444s~ 1 and D2 = 2ir- 6666s -1 , the time of the first returning 
of the wave packet to site 1 becomes, see Fig. 2, 



, R = M_IU 14 . 5 „„, 
The traveling waves of the spin polarization / {li, z ){0) are shown on Fig. 3. 



(32) 
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Figure 3: Propagation of the spin-wave packet along the 201-site open chain. The packet starts at 
boundary site j = 1. 
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4 Intensities of MQ coherences in open chains 
with alternating couplings 

The exact spectrum of Eq.-s ()22l) - (j26j) of the Hamiltonian H (JSJ) provides us with the 
technique for determining of the multi-quantum dynamics in an alternating open chain. 
Again, we take the initial density matrix p(0) ()27|) and calculate how the MQ coher- 
ences develop in the spin system of the alternating spin-^ chain. MQ NMR dynamics of 
the nuclear spins coupled by nearest neighbour dipolar interactions are described by the 
Hamiltonian ^Ql> jllj ■ 

1 

HmQ = ^ X/ Ai,n+1 {In,+In+1,+ + In,— Wl,-} • (33) 
n=l 

The Hamiltonian H ()33|) takes the form of the exactly solvable Hamiltonian H Q (with 
the Larmor frequencies oo n = for all sites) by making use of the unitary transformation 
acting on the even sites, 

Y = exp(-iTrI 2tX ) exp(-iTrI^ x ) ■ ■ ■ exp(-iirI N _ ltX ), (34) 

so that YHmqY + = H( Eq.Q; {w n = 0}). In addition, the transformation Y brings the 
initial density matrix (|27|) to the form 

N 

p(o) = y/ 2 y+ = ^(-ir- 1 / n , 2 , (35) 

n=l 

where we introduce the total polarization I z = J2n=i In,z- The Liouville-von Neumann 
Eq. I)28|) with Hamiltonian H l)33|) and the initial density matrix (|27|) gives the intensities 
G n {t) of n = and n = ±2 orders, just as in the case of the homogeneous chain with the 
conservation condition |13j . 

Go(t) + G 2 (t) + G_2(t) = l, (36) 

and 

, iV N 

G (t) = -J2 cos2 ( X ^ G ±2 (t) = — ^sin 2 (A^). (37) 

n=l n=l 

Fig. 4 demonstrates the development of the 2Q coherence on the alternating N = 201- 
chain with couplings D\ = 2ir ■ 4444 s _1 and D 2 = 2tt ■ 6666 s" 1 , thus 5 = j^- = 1.5. As 
time proceeds, the regular course of the intensities G±2{t) is transformed to the erratic 
temporary behaviour just as in the case of the homogeneous lattice with D\ = D2 ■ F° r 
D± = D2, Eq. (|37|) reproduces exactly the results for the intensities of the 2Q coherence 
ofRef. [HI- 

5 Conclusion 

Sensitivity of the NMR spin polarization dynamics to the spatially periodic short-distance 
inhomogeneity of the lattice has been explored in the previous sections relying on the 
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Figure 4: Time course of the intensities of the 2Q coherence, Eq. Q37|) . in the chain of N = 201 
spins interacting with coupling constants D\ = 2ir • 4444 s _1 and = 2ir ■ 6666 s~ 1 and zero Larmor 
frequencies. Insert shows an early development of 2Q coherence over times t < 1.5 ms. 



exact spectrum of the spin-g XY-Hamiltonian on the alternating open chain with an odd 
number of sites. 

If the spin polarization is prepared initially on the single site, the spin dynamics on the 
non-homogeneous lattice is represented as the propagation of the spin wave packet. The 
velocity of the wave packet on the non-homogeneous lattice is controlled by the minimal 
coupling constant min{Di, Z^}- In the extreme situation D\ <C D2 (D\,D2 7^ 0), the 
spin polarization stays fixed on the initial site. 

As time proceeds, the regular spin propagation along the alternating chain is trans- 
formed to the erratic one. 

The time scale of the regular behaviour of the 2Q intensities G±2(t) is 4 times shorter 
than the time scale of the regular behaviour of the polarization (Ij z )(t)/(Ii z)(0), as it 
can be seen by a comparison of Fig. 4 with Fig. 2. The same effect happens in the case 
of the spin dynamics on the homogeneous chain The effect is caused by two reasons. 
First, the velocity of the propagation of the 2Q coherence is doubled as compared with the 
velocity of the propagation of the spin polarization, as it is obvious by the comparison of 
the dispersion law in Eq. (|3*7j) with the dispersion law in Eq. (|3"U|) . Secondly, the two-spin 
local excitations of the 2Q coherence travel the path iV before the two-spin excitations 
return back on the iV-site chain. The path N should be compared with the path 2N 
traveled by the single local spin excitation before it reappears on the initial site of the 
chain. 
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